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Abstract 

Based on the general form of the master equation for open quantum sys- 
tems the tunneling is considered. Using the path integral technique a simple 
closed form expression for the tunneling rate through a parabolic barrier is 
obtained. The tunneling in the open quantum systems strongly depends on 
the coupling with environment. We found the cases when the dissipation pro- 
hibits tunneling through the barrier but decreases the crossing of the barrier 
for the energies above the barrier. As a particular application, the case of 

decay from the metastable state is considered. 
PACS: 03.65.-w, 05.30.-d, 24.60.-k 
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There has been considerable interest to the quantum tunneling of a particle through an 
energy barrier when the dissipation is presentBHll. Using various models for the description 
of the quantum open system, the opposite dependences of tunneling rate on the dissipation 
have been observed. It is generally thought that tunneling probability decreases in the 
presence of coupling to the environment. Disregarding the stage of averaging over the 
intrinsic degrees of freedom, one can consider the tunneling effect starting right away from 
the general Markovian master equation for the reduced density matrix of the collective 
degree of freedomlllHl] 
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Here, Hq is the Hamiltonian of the collective subsystem and Vj are operators acting in the 
Hilbert space of the subsystem. The second term in (|l|) is responsible for the friction and 
diffusion and supplies the unreversability in the open quantum system. Omitting this term 
we get a standard form for the density matrix evolution equation in the case of closed system. 
The generality of Eq. (|I|) was mathematically proved 

In the one-dimensional case the phase space path integral expression^! for the propagator 
corresponding to ([[]) is written as 

q(t) q'(t) 

G(q, q', t; q , q , 0) = J D[a] J D[a'\ exp (^S[a, a'}), 

qo (o) 9 ; {0 ) 

S[a,a'} = J dT{q(r)p(r) - H ef f(q{T),p(r))} 

dT{q'(T)pi(T)-H: ff (qi(T),p'(T))} 
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with the effective Hamiltonian 

H e ff = H Q - - \ v j\ 2 - 
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Here, the Wigner transform of operators H V^Vj, and Vj are denoted by H , \Vj\ 2 , V* 
and Vj, respectively. For the inverted harmonic oscillator with the Hamiltonian 

and linear environment operators 

V^Atf + Bjq, Vf = A*p + B*q, j = l,2, (4) 
the propagator (0) can be evaluated analytically: 

G(q, q', t; q , q' Q , 0) = 27r fcXh(u;t) exp ^ ex P( iS R/ n ) exp(-Sj/7i), (5) 

where 

- 2cosh(At)[g g - go?'] - 2 sinh(At) [q q' - %q\), 

Si = 8X{u 2_v)smh\u;t) {A{t)[q0 ~ ^ " 4 ^ )[g0 ~ * ~ ? ' ] ~ A( ^~ t)[q ~ ^ 
.4(£) = aexp(2Ai) + 6sinh(2a;t) + ccosh(2cjt) - d, 

= acosh(a;t) sinh(At) + fesinh(u;t) cosh(At), 
■ 2 n r, .2 o\2 



(m^Ja/ - 2A ] - 2XmD pq - D„), 
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- — (m 2 uj 2 D qq + 2XmD pq + D ppj . 
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c = - — (m uj \D qq + XDpp + 2ma; D pq ), 



d=^-{u 2 - \ 2 ){m 2 uj 2 D qq - D pp ). 
Here, the quantum mechanical diffusion coefficients D qq = | £ IA?'I 2 ' D pp = 1 2 l-^il 2 an d 



D gp = -§Re£-4J-B? and the frictional damping rate A = -Im £ A^BjESESB-Ea satisfy the 
following constraints: D w > 0, _D pp > and D pp D qq — D 2 > X 2 h 2 /A which secure the 
non-negativity of the density matrix at any time. The diffusion models, in which these 
constraints are not fulfilled, can be related to the classical or semiclassical considerations 
because they allow the violation of the uncertainty inequality at some 



Using (H), p{t) is determined from p(t = 0) as 



< q\p(t)W >= J d % J dq G(q,q',t;q ,q' ,Q) < q \p(t = 0)\% > . 



(6) 



In order to study the tunneling with the Hamiltonian (Q), we consider a particle in the initial 

state 



*(?) = (2,T(T„(0))- 1/4 exp 



in the left-hand side from the potential barrier. The calculation of (||) with 
yields the Gaussian distribution at time t 
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p(q,t) =< q\p{t)\q >= {2 7 ia qq {t))- 1 / 2 exp( 



-(?-?(*)n, 



(7) 
and (0) 

(8) 



2a qq (t) 

with the first q(t) and second <J qq (t) moments. The equations for these moments are given 
in and below for arbitrary potential. Originally they contain the friction 

in both coordinate \ q and momentum A p so that X p + X q = 2X. The considered particular 
case of A p = \ q = A is generalized for A p ^ \ q by using the canonical transformations^! 
p' = p + pmq, q' = q with the parameter p. Therefore, the expression (H) can be applied to 
the case of A p 7^ X q as well. 

The solutions of equations for the first and second moments in (§) are 



q(t) = e- M ^(0) 
1 



cosh(^/2) + Ap Xq sinh(V>t/2) 



+ — -p(0)sinh(^/2) 
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m(X q — X p ) 2moo'- 
i(X q -X p )C 3 + C 2 ^] sinh(^) 



-\{X q -X p )C 2 + C^} cosh(^) 
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where the following notations are used: 
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and ^ = y (A p — A q ) 2 + 4co> 2 . With these expressions we obtain the same result as in Ref.E3 
at Ap = X q = 0, D pp = D qq = D pq = and a P p(0) = h 2 / (4<r OT (0)) (<x 9 p(0) = 0). For A g = 0, 
D pp = D qq = D pq = and o" P p(0) = h 2 / (Aa qq (Q)) (a qp (0) = 0), our results coincide in 
the underdamped limit with the results of Ref.i where the tunneling was studied with the 
inverted Caldirola-Kanai Hamiltonian. 

The penetration probability at time t is determined by the following expression (g = 
corresponds to the top of the barrier): 

oo 

P(t) = J dq[p(q, t) - p(q, t = 0)]/ J dqp(q, 0), (10) 

-oo 

which is the ratio of change of the probability to be on the right-hand side of the barrier 
in time t over the initial probability of the finding the particle on the entry left-hand side. 
Using Eqs. (p|)-(|l0|), the penetration probability P = P(t — > oo) is easially calculated taking 
the initial variances in accordance with the uncertainty relation. Here, we use a qq (0)a pp (0) = 
h 2 /A and a pq {0) = 0. 

The dependences of the penetration probability through the parabolic barrier on the 
initial energy E of system are presented in Fig. 1 for three sets of the friction coefficients 
Ap and X q . All diffusion coefficients depend only on A. For the sub-barrier energies (E < 
0), the tunneling is larger for X q = X p ^ as comparable to the case without friction 
A P = X q = 0. For E < 0, the dissipation in coordinate X q increases but dissipation in 
momentum A p decreases the barrier penetration. The increase of the tunneling was obtained 
in the microscopic Gisin's modefl for large friction. However, in this model one can not 
distinguish the influence of frictions in coordinate and momentum on the tunneling. Larger 



penetration of the barrier than in the standard coupled- channel calculations is necessary to 
explain the experimental data on the sub-barrier fusion H. It could be that in this case the 
coupling with environment leads to X q ^ that renormalizes the barrier and increases the 
penetrationS. The friction and diffusion reduce the crossing of the barrier for the energies 
above the barrier. For E = and X p = X q , the penetration and reflection probabilities are 
equal to each other with and without dissipation. 

In Fig. 2 we show how the tunneling depends on the diffusion coefficients at different 
values of friction in momentum X p for X q = 0. One can see that only with the diffusion 
coefficient in momentum D pp (D qq = D pq = 0) P decreases with increasing X p . Note that this 
set of diffusion coefficients is not compatible with the quantum mechanical consideration. For 
D qq 7^ 0, the value of P initially decreases with increasing X p up to some "critical" friction 
coefficients and then it starts to grow. This effect becomes more evident at larger D qq 
and D pp (higher temperature). The "critical" friction coefficient decreases with increasing 
temperature. This behaviour of the tunneling probability P as a function of X p can be 
explained in the following way: The tunneling is more crucial to the value of D qq then to the 
value of D pp because cr qq (t) (correspondingly p(q,t) and P) is more sensitive to D qq than to 
D pp ; At large X p the system has a longer time for the tunneling and during this time u qq (t) 
and P(t) strongly increase due to diffusion in coordinate. The increase of tunneling rate 
with temperature is in agreement with RefM 

The probability of finding the particle to the right of the barrier is very sensitive to the 
width <J qq (0) of the initial wave packet localized to the left of the barrier at t = (Fig. 3). 
This effect is weaker with the dissipation. For smaller a qq (0), the value of <J PP (0) becomes 
larger in quantum mechanics and the penetration probability increases due to the larger 
fluctuation energy. In the vicinity of cr qq (0) = Ti/{2mui) the dependence of P on <r qq (0) 
becomes weak and the curve in fig. 3 has a step-like behaviour. 

The calculated time of decay from the metastable state in the potential 

U(q) = aq 2 - (3q 3 (11) 
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is shown in Fig. 4 as a function of X p . These data result from the solution of equations for 
the first and second moments (obtained from eq. ([I])): 



— = -A -+ — 
dt q m 



d j- t = -dU(q)/dq-^d 3 U(q)/dq 3 a qq - X p p, 

— 77^ = —2X q a qq H (jpq + 2D qq , 

dt m 

da, 



JjZ = -2X p a pp - 2d 2 U(q)/dfa pq + 2D pp . 



d U{q)/dq z a qq + — a pp - (X p + X q )a pq + 2D pq . (12) 

These equations are obtained from ([I]) and @ for arbitrary potential U(q). In order to 
calculate P(t) for short times, we can use in the first approximation the formalism elaborated 
for the parabolic barrier. The value of time t\i2 at which P{t\/2) = 0.5 (the value of q 
corresponds to the top of the barrier) may be defined in some sense as the tunneling time@. 
The tunneling time increases monotonically with X p (X q = 0) when D qq = 0. For D qq ^ 0, 
the value of ti/2 initially increases with X p and then it decreases. This means that for large 
X p the dissipation prohibits the decay from metastable state due to diffusion in coordinate. 
The results of calculations above are in agreement with the results obtained in@ using the 
Gisin modelil for the double well potential. 

In conclusion, our calculations show that the dissipative effects on the tunneling pro- 
cess are quite complicated. It is evident that the earlier conclusions that the dissipation 
inhibit tunneling is not correct in the general case. There are examples when the dissipa- 
tion prohibits the penetration through the barrier. Using the general master equation ([!]) 
for describing the open quantum systems, we can transparently show the influence of each 
friction and diffusion coefficient on the tunneling. However, the microscopical calculation 
of these coefficients in the real system remains to be interesting problem. In the consistent 
quantum treatment the tunneling should be calculated with the set of the diffusion coeffi- 
cients where D qq > 0. As was shown, the tunneling is crucial to the value of D qq > 0. If the 
environment operators lead to X q 7^ then the interaction with environment renormalizes 



the potential barrier and influence the tunneling. With the initial Gaussian distribution (^) 
the distribution function remains to be Gaussian at any time. 
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FIGURES 

FIG. 1. Calculated dependence of the penetration probability through the parabolic barrier on 
the initial energy of particle E at temperature kT = 0, hu = 2.0 MeV, q(0) = —2 fm, cr qq (0) = 0.2 
fm 2 , m = 53mo (mo is the mass of nucleon), D qq = hX/(2mu), D pp = Xmhuj/2 and D pq = 0. 
The results for the cases (\ p = X q = 0), (h\ p = HX q = 1 MeV) and (fiX p = 2 MeV, X q = 0) are 
presented by solid, dotted and dashed lines, respectively. 

FIG. 2. Calculated dependence of the tunneling probability P on the friction coefficient 
in momentum X p at temperatures kT = and 3 MeV, fioj = 2.0 MeV, q(0) = —2 fm, 
cr qq (0) = 0.2 fm 2 , E = —5 MeV, m = 53mo, and friction coeffisients in coordinate X q = 
(A = X p /2). The calculations for the cases (kT = 0, D qq = hX/(2muj), D pp = Xmfiui/2 
and D pq = 0), (kT = 0, D qq = 0, D pp = X p mhuj/2 and D pq = 0), (kT = 3 MeV, 
D qq = hX/(2muj)coth(huj/(2kT)), D pp = Xmfioj /2cot\i(fiuj / (2kT)) and D pq = 0) and (kT = 3 
MeV, D qq = 0, D pp = Xpmhuj/2 coth(huj / (2kT)) and D pq = 0) are presented by solid, dashed, 
dotted and dashed-dotted lines, respectively. 

FIG. 3. Calculated dependence of the penetration probability P on the initial variance cr qq (0) 
at kT = 0,Huj = 2.0 MeV, q(0) = -1 fm, p(0) = 0, m = 53m , D qq = hX/(2muj), D pp = Xmhoj/2 
and Dp q = 0. The results obtained with A = A p = X q = and TiX = hX p = %X q = 1 MeV are 
presented by solid and dashed lines, respectively. The value cr qq (0) = h/(2muj) is marked by arrow. 

FIG. 4. Calculated dependence of the decay time from the metastable state in the potential 
(O) on the friction coefficient A p at X q = (A = Xp/2), kT = 0, p(0) = 0, a qq (0) = 0.2 fm 2 



and m = 53mo- The depth of potential pocket with the minimum at g(0) = —1.08 fm is 5 MeV 
(a = —2.57 MeV fm -2 and = 1.59 MeV fm -3 ). The top of the barrier corresponds E = MeV 
at q = fm. The calculations for the cases (D qq = hX/ (2muj), D pp = Xmfiuj/2 and D pq = 0) and 
(D qq = 0, D pp = X p mhu)/2 and D pq = 0) are presented by solid and dotted lines, respectively. 
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